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We present an example of a partial dynamical symmetry (PDS) in an interacting 
fermion system and demonstrate the close relationship of the associated Hamiltonians 
with a realistic quadrupole-quadrupole interaction, thus shedding new light on this 
important interaction. Specifically, in the framework of the symplectic shell model 
of nuclei, we prove the existence of a family of fermionic Hamiltonians with partial 
SU(3) symmetry. We outline the construction process for the PDS eigenstates with 
good symmetry and give analytic expressions for the energies of these states and 
E2 transition strengths between them. Characteristics of both pure and mixed- 
J> ' symmetry PDS eigenstates are discussed and the resulting spectra and transition 

strengths are compared to those of real nuclei. The PDS concept is shown to be 
<^ ■ relevant to the description of prolate, oblate, as well as triaxially deformed nuclei. 

Similarities and differences between the fermion case and the previously established 
partial SU(3) symmetry in the Interacting Boson Model are considered. 

o 

^ ; PACS numbers: 21.60.Fw, 21.10.-k, 27.20.+n, 27.30.+t 
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I. INTRODUCTION 



- 1—^ 1 

^ ■ Symmetries play an important role in physics. Constants of motion associated with a 

symmetry govern the integrability of a given classical system, and at the quantum level 
symmetries provide labels for the classification of states, determine selection rules, and 
simplify the relevant Hamiltonian matrices. Algebraic, symmetry-based, theories have been 
firmly established as an elegant and practical approach to a variety of physical systems 
(see, for example, Refs. ||, §, |, g §, g g §, |, 0, 01], 0). These theories offer the 



greatest simplifications when the interaction under consideration is symmetry-preserving in 
the selected state labeling scheme, that is, when the Hamiltonian either commutes with all 
the generators of a particular group ('exact symmetry') or when it is written in terms of and 
commutes with the Casimir operators of a chain of nested groups ('dynamical symmetry'). 
In both cases basis states belonging to inequivalent irreducible representations of the relevant 
groups do not mix, the Hamiltonian matrix has block structure, and all properties of the 
system can be expressed in closed form. An exact or dynamical symmetry not only facilitates 
the numerical treatment of the Hamiltonian, but also its interpretation and thus provides 
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considerable insight into the physics of a given system. 

Naturally, the application of exact or dynamical symmetries to realistic situations has its 
limitations: Usually the assumed symmetry is only approximately fulfilled, and imposing 
certain symmetry requirements on the Hamiltonian might result in constraints which are too 
severe and incompatible with experimentally observed features of the system. The standard 
approach in such situations is to break the symmetry. In cases where a symmetry-breaking 
Hamiltonian is involved, it is possible to decompose the offending terms into basic parts 
("irreducible tensor operators") which exhibit specific transformation properties. Provided 
the appropriate group coupling coefficients and the matrix elements of some elementary 
tensor operators are available, matrix elements of operators which connect inequivalent 
irreducible representations can be determined and the exact eigenvalues and eigenstates can 
then be obtained (at least in principle). While group theoretical considerations still play 
an important role in evaluating the coupling coefficients and matrix elements for such a 
calculation and in truncating model spaces which have become too large for a complete 
numeric treatment, the basic simplicity of the symmetry-based approach is lost. 

Alternatively, one might consider some intermediate structure, which allows for symmetry 
breaking, but preserves the advantages of a dynamical symmetry for a part of the system. 
Partial dynamical symmetry (PDS) |L3| provides such a structure. It corresponds to a 
particular symmetry breaking for which the Hamiltonian is not invariant under the symmetry 
group and hence various irreducible representions (irreps) are mixed in its eigenstates, yet 
it possesses a subset of 'special' solvable states which respect the symmetry. The notion of 
partial dynamical symmetry generalizes the concepts of exact and dynamical symmmetries. 
In making the transition from an exact to a dynamical symmetry, states which are degenerate 
in the former scheme are split but not mixed in the latter, and the block structure of the 
Hamiltonian is retained. Proceeding further to partial symmetry, some blocks or selected 
states in a block remain pure, while other states mix and lose the symmetry character. 

Other generalizations of the idea of dynamical symmetry are possible. Van Isacker |T3 



for example, suggested to break the dynamical symmetry associated with an intermediate 
group G2 in a subchain G\ D G2 D G3 for all states of the system, while preserving the 
remaining (dynamical) symmetries. The resulting Hamiltonian is in general not analytically 
solvable, but its eigenstates can still be (partly) classified by quantum labels associated with 
the groups G\ and G 3 . An approximate-symmetry scheme called quasi-dynamical symmetry 
was discussed by Bahri and Rowe [jnj. They considered strong but coherent mixing of the 
irreducible representations associated with a given dynamical symmetry. Both methods of 
extending the concept of dynamical symmetry differ from the notion of partial dynamical 
symmetry introduced above since, unlike in the partial-symmetry case, the eigenvalues of 
the Hamiltonians cannot be obtained analytically, not even for a part of the system. 

The partial-symmetry scheme was introduced first in bosonic systems, where it was ap- 



plied to the spectroscopy of deformed nuclei. In Ref. [16|, a Hamiltonian with partial SU(3) 
symmetry was constructed in the framework of the Interacting Boson Model (IBM) of nuclei 
||, and the calculated spectrum and E2 rates of 168 Er were compared to experimental re- 
sults. The PDS Hamiltonian was found to reproduce the experimentally observed feature of 
non-degenerate rotational 7 and (3 bands ('K-band splitting') and to possess several bands of 
solvable states, whereas previous attempts to describe the 168 Er data had involved Hamilto- 
nians with SU(3) dynamical symmetry, which can only yield 7 and j3 bands with degenerate 
angular momentum states, or had achieved agreement with the data by completely break- 
ing the SU(3) symmetry. Employing the same Hamiltonian, Sinai and Leviatan |T7, I? 
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investigated the structure of the lowest collective K=0 + excitation in deformed rare-earth 
nuclei. Implications of the partial dynamical symmetry for the mixing behavior of this col- 
lective band were discussed and compared to broken-SU(3) predictions. In another study 



Ref. JOl, in the context of the IBM-2, the proton-neutron version of the Interacting Boson 



Model U 20], Talmi was able to explain simple regularities in spectra of the Majorana op- 



erator as an example of partial dynamical symmetry. More recently, the relevance of partial 
F-spin symmetry was studied in the framework of the IBM-2. It has long been known that 
F spin, the SU(2) quantum number associated with the two-valued proton-neutron degree 
of freedom of the IBM-2, cannot be conserved in nuclear spectroscopy. However, Leviatan 
and Ginocchio demonstrated that empirical energy systematics in the deformed Dy-Os 
region can be reproduced under the assumption of partial F-spin symmetry. Moreover, the 
associated partial symmetry Hamiltonians point to the existence of F-spin multiplets of scis- 
sors states, with a moment of inertia equal to that of the ground band. These predictions 
were tested against recent analyses of Ml transition strengths. 

The subject of partial symmetries and supersymmetry in nuclear physics was considered 
by Jolos and von Brentano in the context of the Interacting Boson-Fermion Model |2lJ and 
the particle-rotor model Q . 



Partial symmetries can be associated with continuous as well as discrete groups. The 



dynamical groups employed in the IBM, e.g., are continuous. In Ref. p4| , an example of 
a partial symmetry which involves point groups was presented in the context of molecular 
physics. Ping and Chen used a model of N coupled anharmonic oscillators to describe the 
molecule XYq. The partial symmetry of the Hamiltonian allowed them to derive analytic 
expressions for the energies of a set of unique levels and to discuss the structure of the asso- 
ciated eigenstates. Furthermore, the numerical calculations required to obtain the energies 
of the remaining (non-unique) levels were greatly simplified since the Hamiltonian could be 
diagonalized in a much smaller space. 

Partial symmetries have relevance not only for discrete spectroscopy but also for the study 
of stochastic properties of dynamical systems. A generic classical or quantum-mechanical 
Hamiltonian exhibits mixed dynamics: areas of regular motion and chaotic regions coexist 
in phase space, and even when a system seems to be fully chaotic, regular states may 
exist. Whelan et al. |2^j used Hamiltonians with partial dynamical symmetries to investigate 
quantum-mechanical systems which are partly regular and partly chaotic. In the context 
of the Interacting Boson Model, it was demonstrated that partial symmetries impose a 
particular phase-space structure which leads to a suppression of chaos in mixed systems. 



Canetta and Maino |26| carried out a quantum-statistical analysis of regular and chaotic 
dynamic behavior in the IBM-2. Varying the Hamiltonian parameters, they observed a 
nearly regular region in parameter space — far away from dynamical symmetry limits - 
which they linked to the existence of a partial dynamical symmetry. Since Hamiltonians 
with partial symmetries are not completely integrable and may exhibit stochastic behavior, 
they are an ideal tool for studying mixed systems with coexisting regularity and chaos. 

Partial symmetries are not confined to bosonic systems. In Ref. ||27|| , an example of a 
partial symmetry in an interacting fermion system was presented. A family of Hamiltonians 
with partial SU(3) symmetry was introduced in the framework of the symplectic shell model 



of nuclei [^8| . The Hamiltonians were shown to be closely related to the deformation-inducing 
quadrupole-quadrupole interaction and to possess both mixed-symmetry and solvable pure- 
SU(3) rotational bands. For the example of the (prolate) deformed light nucleus 20 Ne, it 
was demonstrated that various features of the quadrupole-quadrupole interaction can be 
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reproduced with a particular parametrization of the PDS Hamiltonians. In that work, the 
partial dynamical symmetry was identified directly at the fermion level. It is also possible 
to start with a bosonic PDS Hamiltonian and map the bosonic generators into fermionic 
generators of the same algebra. This approach was taken by Mamistvalov |2{J , who studied 
partial symmetry in a schematic SU(2) xSU(2)-type Lipkin model. Very recently, partially 
solvable shell-model Hamiltonians with seniority-conserving interactions were investigated 



by Rowe and Rosensteel p0|. 



It is the purpose of this work to investigate the fermionic PDS Hamiltonians presented 
in Ref. |27| in more detail. Specifically, the construction process for the pure eigenstates 



is outlined and analytic expressions for the energies of pure states and the strengths of 
E2 transitions between these states are given. Properties of the special solvable states are 
discussed and an application to the oblate deformed light nucleus 12 C and the prolate nucleus 
20 Ne are presented. Moreover, an application to 24 Mg demonstrates the relevance of the PDS 
concept for well-deformed, triaxial nuclei. 



In the next section, the symplectic shell model (SSM) is reviewed. In Section HI, a family 
of symplectic Hamiltonians with partial SU(3) symmetry is introduced and their relation to 
the quadrupole-quadrupole interaction is established. Properties of the special eigenstates 
of the PDS Hamiltonians are discussed in Section [TV| , and applications to realistic nuclear 
systems are presented in Section [V|. In Section VI, the fermionic PDS Hamiltonians are 



compared to the earlier introduced bosonic PDS Hamiltonians [0, [17], [Tj|, and Section |VT l 
summarizes our work. Appendix A contains further relevant material regarding SU(3) cou- 
pling coefficients and reduced matrix elements and Appendix B presents expressions for 
matrix elements of operators employed in the calculations. 

II. THE SYMPLECTIC SHELL MODEL 

The symplectic shell model (SSM) is an algebraic, fermionic, shell-model scheme which 
includes multiple 2fru one particle-one hole excitations. It includes all essential observables 
for a description of nuclear monopole and quadrupole collective vibrations as well as for rigid 
and irrotational flow rotations. Since the model allows for intershell excitations and since 
its observables are expressible in microscopic shell-model terms, it provides a multi-shell 



realization of the nuclear shell model [28 



A. Symplectic Generators 

The symmetry algebra of the symplectic scheme is spanned by one-body operators which 
are bilinear products in the (relative) position (x S j, % = 1,2,3, s = 1, ...,A— 1) and 
momentum (p S i) observables: 



Qij ^ ^ X s iX s j 

s 

K <3 = ^VsiVsj (1) 

s 
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Sij ^ ^ tp^siPsj Psi-E sj) j 

s 

where A — 1 is the number of Jacobi 'particles' remaining after removal of the center-of- 
mass contribution. Together the operators generate the 21-dimensional symplectic algebra 
sp(6,R), that is, the Lie algebra of linear transformations which preserve a skew-symmetric 
bilinear form on a six-dimensional real vector space. It is the smallest Lie algebra that 
contains both the quadrupole moments and the many-nucleon kinetic energy, and it has 
several physically relevant subalgebras. These include gcm(3) and its subalgebra [i? 5 ]so(3), 
associated with the Geometric Collective Model and its rotational limit, respectively, the 
algebra gl(3, R) of the general linear motion group, as well as su(3) and its subalgebra 
so(3), associated with the Elliott model and the rotation group, respectively. The sp(6,R) 
algebra furthermore includes the canonical subalgebras sp(2,R) and sp(4,R), which have 
been studied by Arickx et al. f3lfl , and by Peterson and Hecht |32| , respectively, as possible 



approximations to the full three-dimensional symplectic model. 

For many purposes, it is advantageous to express the symplectic generators in terms 
of harmonic oscillator boson creation and annihilation operators b\ { = (x S i — ip S i)/\/2 and 
b s = (x S i+ip S i) I a/2. The symplectic generators may then be expressed as one-body operators 
which are quadratic in the oscillator bosons |33 : 



s 

= \ Y. h ^ ( 2 ) 

s 



Alternatively, one may use the spherial components of the oscilator bosons, b 



s,l,±l 



T^g(&Ii ± ib\ 2 ), 6$°} = b\ 3 , and bf£ ±1 = T^(b s i ± ib s2 ), bf£ = b s3 , to write the gen- 
erators as SU(3) tensor operators pi R3]: 



v(H) _ ./or rz,t(io) v I(oi) ,(11) 



<5£f = V2E,{i.r ) x fT},™ (.1 = 1,2) 

(3) 

II, ., J,' I'-"' ,■; n ■)! 

■J2 



A? = ^E,{d (10, xd (10) }i;i 0) « = o,2) 



The notation indicates that the operator T possesses good SU(3) [superscript (A,/i)] 

and SO (3) [subscript Im] tensorial properties. Since b\b\ adds two quanta to particle s, 
thereby moving it up across two major oscillator shells, A^' creates a 2ftw excitation in 
the system. Analogously, which is related to Aj^ by Hermitean conjugation, = 

(_1 ) l ~ m (Af^ n )\ annihilates a 2ftw excitation. The act only within a major harmonic 
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oscillator shell. They generate the group SU(3) of the well-known Elliott model [36]: 
V3C^ = Q E 2m = ^^(r s 2 F 2m (f s ) + p% m (p s )) (m = 0, ±1, ±2) , 



(4) 



A(n) 



L q (g = 0,±l), 



where Q^m denotes the symmetrized quadrupole operator of Elliott, which does not couple 
different major shells, and L q is the orbital angular momentum operator. The harmonic 
oscillator Hamiltonian, H = Yli=i ls a SU(3) scalar and generates U(l) in U(3) = 
SU(3) x U(l). 

Alternatively, one can realize the symplectic generators in terms of fermionic creation 
and annihilation operators [31|: 



,(11)5=0 



(20)5=0 



(5) 



B 



(02) 
Im 



E v \/Mv + 1){V + 2)(7? + 3)(r7 + 4){at x a v+2 }^Z°o + OT(A) 



Wllere a llml/2a K*ml/2 



l)' 7+ ' +m+1 / 2+0 'a r/ /_ m i/2 -a) is a single-particle creation (annihi- 
lation) operator, which produces (destroys) a fermion with angular momentum I, projection 
m and spin 1/2, projection a in the 77-th major oscillator shell. The sums run over all 
shells, and the coupling to total spin 5 = with projection S = reflects the fact that the 
generators act on spatial degrees of freedom only. The operators O cm (A) remove the spurios 
center-of-mass content from the generators. Details regarding the fermionic realization of 



Sp(6,R) can be found in Ref. |33] 



B. Symplectic Basis States 



A basis for the symplectic model is generated by applying symmetrically coupled products 
of the 2hcu raising operator A- 2 ^ with itself to the usual Ohio many-particle shell-model 
states. Each Ofuv starting configuration is characterized by the distribution of oscillator 
quanta into the three cartesian directions, {01, a 2 , 03}, where <J\ > a 2 > cr 3 . Here <7j denotes 
the eigenvalue of the U(3) weight operator Ca = ^ s (&sj& S i + 1/2), which essentially counts 
the number of oscillator bosons in the i-th direction of the system. Since s — 1, 2, . . . , A — 1, 
it follows that the cij are half-integer numbers for even-A and integers for odd-A nuclei. 
Equivalently, one may employ quantum numbers N a (X a , /i CT ), where A CT = <J\ — a 2 , \i a = 
o~2 — 03 are the Elliott SU(3) labels, and N a = o\ + cr 2 + cr 3 is the eigenvalue of the harmonic 
oscillator Hamiltonian, ifo = C11 + C22 + C33, which takes the minimum value consistent with 
the Pauli Exclusion Principle. Each such set of U(3) quantum numbers uniquely determines 
an irreducible representation (irrep) of the symplectic group, since it characterizes a Sp(6,R) 
lowest weight state. Any component of the symplectic lowering operator (and of 

with i < j) annihilates such a lowest weight state. 
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In contrast, application of the symplectic generator allows one to successively build 
a basis for the Sp(6,R) irrep under consideration: The product of N/2, N = 0,2,4,..., 
raising operators A*- 20 ) is multiplicity-free and generates Nfrw excitations for each starting 
configuration N a (\ a , fx a ). Each such product operator p 7V ( A ™'^™) can be labeled according to 
its U(3) content, {ni,ri2,n 3 } or N(\ n ,[i n ), where (A n ,/z n ) ranges over the set 



{{ni — n 2 , n 2 — n 3 )\rii > n 2 > n 3 > 0] N = rii + n 2 + n 3 ; ri\, n 2 , n 3 even integers} (6) 



The raising polynomial 'p N ( x "->Vn) is then coupled with \N a (X a , /x CT )) to good SU(3) symmetry 
/o(Ao), Hu,), with p denoting the multiplicity in the coupling (A n ,/i n ) (g> (A CT ,/i CT ). The quanta 
distribution in the resulting state is given by {uj\, u 2 , u 3 }, with N^ = N a + N = oj\ + u 2 + u 3 , 
^1 > ^2 > k>3, and Xuj = lo\ —u> 2 , /j^ = uj 2 — uj 3 . The states of the Sp(6,R) D SU(3) basis are 
thus labeled by three types of U(3) quantum numbers: T a = {a\, 0%, a 3 } = N a (X a , fi a ), the 
symplectic bandhead or Sp(6,R) lowest weight U(3) symmetry, which specifies the Sp(6,R) 
irreducible representation; T n = {rii,n 2 ,n 3 } = iV(A n // n ), the U(3) symmetry of the raising 
polynomial; and T w = {ui,u 2 ,u 3 } = N u (\ u , (j, u ), the U(3) symmetry of the coupled prod- 
uct. A given symplectic representation space N a {\ ai jji a ) is infinite dimensional, since N/2, 
the number of oscillator excitations, can take any non-negative integer value. In practical 
applications, one must therefore either truncate the symplectic Hilbert space, or restrict 
oneself to interactions and observables for which the matrix elements depend solely on the 
symplectic irrep and can be calculated analytically. The basis state construction is schemat- 
ically illustrated in Fig. [I] for a typical Elliott starting state with (A CT , = (0, /x). A similar 
figure for (\ a ,fj, a ) = (A, 0) is given in Ref. P7J. 



To complete the basis state labeling, additional quantum numbers a are required. This 
can be accomplished by reducing Sp(6,R) D SU(3) states with respect to the subgroup U(l) 
x SU(2) of SU(3) and assigning labels a = eAM\ 1 . This SU(2) subgroup, however, is 
not the physical orbital angular momentum subgroup SO(3) of SU(3). States with good 
angular momentum values can be obtained from the SU(3) D U(l) x SU(2) (canonical) 



basis by projection |36| , |37| . The associated quantum numbers are a = kLM, where k is 
a multiplicity index, which enumerates multiple occurrences of a particular L value in the 
SU(3) irrep (A,/i) from 1 to K^ ax (A, fi), 

kT x {\ /i) = [(A + 11 + 2 - L)/2] - [(A + 1 - L)/2\ -[(// + 1 - L)/2\ , (7) 

where [. . .] is the greatest non-negative integer function [|38|| . The K™ ax (\,fi) occurrences of 
L can be distinguished in a variety of ways. The physically most significant scheme is that of 
Elliott 1 5^| , in which case the projection of L along the body- fixed 3-axis, denoted K, is used 



to sort the L-values into the familiar K-bands of the rotational model. Unfortunately, states 
defined in this manner are not orthonormal with respect to the multiplicity quantum number 
K. To avoid the resulting complications, such as working with non-hermitean matrices, the 
Elliott basis is usually orthonormalized using a Gram-Schmidt process. Vergados |3H[| , for 
example, gives a prescription to construct orthogonal basis states in a systematic manner for 
all (A,/i), such that the physical interpretation of K as a band label can be approximately 



1 Here e, the eigenvalue of Q20, gives the U(l) content and the SU(2) irrep is characterized by A with 
projection M\. 
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maintained 2 . In the present work, we employ the orthonormal basis of Vergados. For 
simplicity, however, we use the running index to distinguish multiple 

occurrences of L in a given SU(3) irrep (A,/i) and list the corresponding Vergados labels 
where appropriate. The dynamical symmetry chain and the associated quantum labels of 



the above scheme are given by [28 
Sp(6,R) 



D U(3) D SO(3) D 50(2) 

N a (X a , fia) N(X n ,II n )p N^XuifJLu) k L M 



(8) 



When applying the formalism to realistic nuclei, we assign rotational band labels according 
to the calculated B(E2) rates. 

The quadratic Casimir operators of SU(3) and Sp(6,R), 



( - / 2 ' u 2 



£,(11) _ £,(11) 



CsU3 — 

C,s P 6 — —2A y ^'B^' — 2A^ U ' ■ B^' + C S U3 + 77-^0 — 4if 
have the following eigenvalues in the dynamical symmetry basis: 
(CW[(A,/i)] = 2(A 2 + /i 2 + A/i + 3A + 3/x)/3 



+ C S U3 + -^Hq 



(9) 
(10) 

(11) 



(Csp^lNAK,^)] = 2(A 2 + X a fi a + 3X a + 3fi a )/3 + N 2 /3-AN a . (12) 

The collection of all Ohu configurations provides a complete Hilbert space for the Elliott 
SU(3) submodel of the SSM and is referred to as the Ofku horizontal shell-model space. 
The set of states built on a given U(3) irrep N a (X a fi a ) is called the vertical extension of 
that irrep. Each vertical extension can be partitioned into horizontal slices with the states 
within the iV/2-th slice representable as a homogeneous polynomial of degree N/2 in the 
tensors acting on the parent OJtlu configuration (see also Fig. ||]). Interactions can 
thus be classified according to their effect on this structure; pairing, for example, causes 
horizontal mixing, both within each 'cone' (symplectic irrep) and between different cones, 
while the quadrupole-quadrupole interaction induces horizontal and vertical mixing, but 
does not connect different cones. 



C. Symplectic Hamiltonians 

A primary goal of the symplectic shell model is to achieve a microscopic description 
of deformed nuclei. These nuclei exhibit collective behavior, that is, modes of excitation 



2 Vergados projects from a state with e — e m in = — A — 2/i, A = A/2, M\ = A/2 for A > fx and 
e = e max = 2A + fi, A = /ii/2, Ma = —fi/2 for A < fi and employs the 'Elliott rule' to determine the 
possible K values, K = min(A, /i), min(A, /i) - 2, . . . , 1 or 0, and angular momenta, L = K, K+l, K+2, . . ., 
K+max(A, /i) for K ^ and L=max(A, /i), max(A,/i)-2, . . ., 1 or for K=0. It is also possible to project 



from e = e ma x, A = /«/2, and Ma = +m/2 or s = s m i n , A = A/2, M\ = —A/2. Draayer et al. |37|. |40 
discuss the different projection possibilities and give rules analogous to the Elliott rule for determining 
the K and L content of a given SU(3) irrep (A,//). 
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in which an appreciable fraction of the nucleons in the system participate in a coherent 
manner, as, for example, is the case for rotations. An appropriate Hamiltonian for describing 
rotational phenomena within the symplectic model consists of the harmonic oscillator, which 
provides the background shell structure, the quadrupole-quadrupole interaction, Q 2 -Q2, and 
a residual interaction that should include, for example, single-particle spin-orbit and orbit- 
orbit terms, as well as pairing and other interactions. However, most applications of the 
theory are much less ambitions than this, restricting the interaction to terms that can be 



expressed solely in terms of generators of the symplectic algebra P3], IT, 031, Interactions 
of the latter form do not mix different symplectic irreps and therefore the Hamiltonian 
matrix for such interactions becomes block-diagonal. Indeed, in most practical applications 
the Hilbert space of the system is truncated to one single symplectic representation. This is 
accomplished by selecting the vertical slice (symplectic irrep) constructed from the leading 
starting irrep of the Ofau space. The leading irrep is defined to be the U(3) representation, 
N a (X a , fi a ), from the lowest layer with the most symmetric spatial permutation symmetry 
consistent with the Pauli principle, and the maximal possible SU(3) Casimir operator value, 
(Csu3) [(A CT , Ha)}- For 12 C, for instance, the leading irrep is given by N a (X a ,fi a ) = 24.5 
(0,4), which corresponds to the symplectic weights a 1 =<T 2 =9.5, <t 3 =5.5; for 20 Ne, one finds 
N a (\ a ,n a ) = 48.5 (8,0), since ^=21.5, a 2 =a 3 =13.5 p§; and 24 Mg has N a (A CT ,/v) = 62.5 
(8,4), that is, <Ti=27.5 02=19.5, o"3=15.5 J42|. The single-symplectic irrep approximation is 



a sensible choice for nuclear systems which have a dominant quadrupole-quadrupole force, 
since this interaction does not mix symplectic representations and favors states with large 
(C S U3)[(K,^a)} values. 

A typical Hamiltonian for a calculation in a space truncated in the manner described 
above, is given by a harmonic oscillator term, H Q , plus a collective potential, and a residual 
interaction: 

H = H + Vcoii + V res . (13) 

We choose the collective potential to be a simple quadratic, rotationally invariant, function 
of the microscopic quadrupole moment 3 , Q 2 m = \J ^ XL r s^2m(^s), namely 

Vcoii = -xQ 2 ■ Q 2 • (14) 

The quadrupole-quadrupole interaction is a standard ingredient in models that aim at re- 
producing rotational spectra and nuclear deformations. It emerges (apart from a constant) 
as a leading contribution in the multipole expansion of a general two-body force. It mixes 
states from different oscillator shells, since the quadrupole operator has non-vanishing ma- 
trix elements between shells differing by zero or two oscillator quanta. A major strength 
of the symplectic model is its ability to fully accommodate the action of the quadrupole 
operator, which can be written in terms of symplectic generators: 

Qi m = V3(C^+AS + B^). (15) 

As a result, the model is able to reproduce intra-band and inter-band E2 transition strengths 
between low-lying, as well as giant resonance, states without introducing proton and neutron 
effective charges. 



3 Higher order rotational scalars in Q2 can be included in V co u in order to accomodate more complicated 
potential forms, e.g. a cubic term introduces a 7-dependence into the potential. 
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The effective residual interaction, V res , is included to replace non- collective components 
of a more realistic Hamiltonian and the neglected effects of couplings to other Sp(6,R) 
representations. As in previous works, we choose V res to be a rotationally invariant function 
of the SU(3) generators. For prolate and oblate nuclei we use: 

V res = d 2 L 2 + dj} , (16) 

where L denotes the angular momentum operator, Eq. ([|). This allows us to reproduce 
the energy splittings between states of a rotational band. For triaxial nuclei, such as 24 Mg, 
it becomes necessary to include further terms, in order to reproduce the experimentally 
observed 'K-band splitting', the energy differences found between states with the same 
angular momentum but different K-band assignments. This can be achieved by including 
'SU(3) D SO(3) integrity basis' operators X 3 = (L x Q E ) {1) ■ L and X 4 = (L x Q E ) W ■ (L x 
Q E )(i) in the residual interaction JPJ: 



V' res = C3X3 + c 4 X 4 + d 2 L 2 + d 4 L 4 . (17) 

The evaluation procedure for the matrix elements of the symplectic generators A^ 20 \ 
B^ 2 \ and C^ ll \ and of the integrity basis operators X 3 and X 4 is discussed in Appendix B. 



III. PDS HAMILTON! ANS AND THE QUADRUPOLE-QUADRUPOLE 

INTERACTION 

In this section we introduce a family of fermionic Hamiltonians with partial dynamical 
symmetry. Motivated by the fact that a realistic quadrupole-quadrupole interaction breaks 
SU(3) symmetry within a given major oscillator shell, we define a family of Hamiltonians, 
H(f3o,/3 2 ), which allows us to study the features of the symmetry-breaking terms in some 
detail. The new Hamiltonians do not couple different oscillator shells and, for a particu- 
lar choice of the parameters (3$ and (3 2 , reduce to a form which is closely related to the 
quadrupole-quadrupole interaction restricted to a shell. We prove that this family of Hamil- 
tonians exhibits partial SU(3) symmetry and give rules for determining the 'special' irreps 
and the associated pure eigenstates. 

In the symplectic shell model, the quadrupole-quadrupole interaction can be expressed 
in terms of symplectic generators |45||: 



Q 2 -Q 2 = 3((7 2 + A 2 + B 2 ) ■ {C 2 + A 2 + B 2 ) . (18) 



Employing the commutation relations B 2 ■ A 2 — A 2 ■ B 2 = yi^o and B 2 -C 2 — C 2 - B 2 = ^g-B, 
given in Ref. [45|, this can be rewritten as: 



Q 2 -Q 2 = 3C 2 ■ C 2 + 6A 2 ■ B 2 + 10# 



+ 



(&C 2 ■ B 2 + 10V6B + 3B 2 ■ 4) + H.c] (19) 



where 3C 2 ■ C 2 = Q E ■ Q 2 and H.c. denotes the Hermitian conjugate of the expression in 
parentheses. The first three terms in the expansion, Eq. (0), act solely within a major 
harmonic oscillator shell, while the second line connects states differing in energy by ±2fkj 
and ±4hu). It is primarily the presence of the mufti- /la; correlations that differentiates the 
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SSM from the Elliott SU(3) model. The symplectic model allows for coherent multi-shell 
admixtures in its wave functions and thus achieves the experimentally observed nuclear 
deformation and absolute B(E2) rates. In contrast, the Elliott model requires effective 
charges, since it employs the algebraic (or Elliott) quadrupole-quadrupole interaction, 

Qf ■ Qi = 6C SU3 - 3L 2 , (20) 

which does not connect different oscillator shells. 

Although matrix elements of Q2 and are identical within a harmonic oscillator shell, 
the corresponding quadrupole-quadrupole interactions exhibit differences here as well: The 
62 ■ C2 and Hq terms in the expansion, Eq. (|P9|), are diagonal in the dynamical symmetry 
basis, Eq. (H), whereas A 2 ■ B 2 contains contributions which mix different SU(3) irreps. This 
follows from the relations: 

(20)A(02)_ 1 rl A, (00) J 5 (A^h\( 22 ) 



A B = Aq 'Bo = ^={A x B}r ~ f^{A x B} ' 

f(20) _ A(02) _ 5 j - (00) /5 « - (22) 



a 2 ■ b 2 s ■ = x s >r + v x 5 >o , (21) 

where 

{i x Eji 00 ) = ^(Csus + \h% - AH - C Sp6 ) . (22) 

The term {A x -Bj-Q 00 '* is a SU(3) scalar, but {A x -BIq 22 -* breaks SU(3) symmetry. Within a 
major oscillator shell, it is mainly this symmetry-breaking term that distinguishes the action 
of Q2 ■ Q2 from the effect of the Elliott interaction, Q§ ■ Qf , which respects the symmetry. 

To explore this latter aspect in more detail, we rewrite the collective quadrupole- 
quadrupole interaction as follows: 

Q2 ■ Q2 = QCsu3 — 3Cs p6 + Hq — 2H — 3L 2 — 6A B 

+{terms coupling different h.o. shells} . (23) 

The quadratic Casimir invariants of SU(3), Csu3, and of Sp(6,R), Csp6, and their eigenvalues, 
are given in Eqs. (|9D-(|T2"D. In order to focus on the action of Q2 ■ Q2 within a harmonic 
oscillator shell, we introduce the following family of rotationally invariant Hamiltonians: 



H(p ,p 2 ) =PoAoB + P2A2-B 2 

- QCspe + 3^ 2 - 36^ ) + (A) - fo)A B . (24) 



U 

For p = j3 2 , one recovers a Sp(6,R)DSU(3) dynamical symmetry Hamiltonian: H(/3 ,j3 2 = 
Po) contains only SU(3)-scalars, that is, it does not mix different SU(3) irreps. Furthermore, 
all eigenstates at a given Nhu excitation which belong to the same symplectic and SU(3) 
irreps are degenerate. Additional SO (3) rotational terms, such as L 2 and L 4 split the 
degeneracies, but do not change the wave functions. For the special choice Pq = 12, p 2 = 
18, one finds that H(p = 12, /3 2 = 18) is closely related to the quadrupole-quadrupole 
interaction: 

Q 2 -Q 2 = H(p = 12, p 2 = 18) + canst - 3L 2 

+{terms coupling different h.o. shells} , (25) 
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where the value of const = 6Cs P 6 — 2H Q 2 + 34:H is fixed for a given symplectic irrep N a (X a , fi a ) 
and Nfrw excitation. Although H(fa, fa) does not couple different harmonic oscillator shells, 

A A (22) 

it contains the SU(3)-symmetry breaking term {A x B} and is therefore expected to 
exhibit in-shell behavior similar to that of Q 2 • Qi- 

From Eq. (|2l|) it follows that H(/3q,/32) is generally not SU(3) invariant. We will now 
show that H(fa,fa) exhibits partial SU(3) symmetry. Specifically, we claim that among the 
eigenstates of H(fa, fa), there exists a subset of solvable pure-SU(3) states, the SU(3)dSO(3) 
classification of which depends on both the Elliott labels (A CT , u a ) of the starting state and the 
symplectic excitation N. In general, we find that all L-states in the starting configuration 
(N = 0) are solvable with good SU(3) symmetry (\ a ,u a ). For excited configurations, with 
N > (N even), we distinguish two possible cases: 

(a) Ao- > n a : the pure states belong to (A, /x) = (A CT — N, u a + N) at the Nftw level and have 
L = He+N^e+N+l, . . .,\ a -N+l with N = 2,4, . . . subject to 2N < (X a -/i (T + l). 

(b) A,j < u a : the special states belong to (A, fi) = (A CT + N, u a ) at the Nhw level and have 
L = A CT + N, X a + N + 1, . . . , K + N + n a with iV = 2, 4, — 

To prove the claim, it is sufficient to show that Bo annihilates the states in question 
(since H((3q = fa) is diagonal in the dynamical symmetry basis). For N = this follows 
immediately from the fact that the Ohut starting configuration is a Sp(6,R) lowest weight 
which, by definition, is annihilated by the lowering operators of the Sp(6,R) algebra. The 
latter include the components of the generator B^ 02 \ 

For iV > 0, we have to consider the action of B in more detail. Let {01, 02, 03} be the 
quanta distribution for a Ofuv state with A CT > \i a . An excited state with SU(3) character 
(A,/i) = (A CT — N, ii a + N) must have the quanta distribution {a\,a 2 + N,a 3 }. Acting 
with the rotational invariant B$ on such a state does not affect the angular momentum, 
but removes two quanta from the 2-direction, giving a (N — 2)fkv state with (A',//) = 
(A CT — N + 2, fi a + N — 2). Note that the symplectic generator B cannot remove quanta 
from the other two directions of this particular state, since this would yield a state which 
has fewer oscillator quanta in the 1- or 3-direction than the starting (Ohuu) configuration, 
i.e. the resulting state would belong to a different symplectic irrep. Comparing the number 
of occurrences of a given angular momentum value L in (A,//) at Nftw and (A',//) at (N — 
2) foui, one finds the following: As long as A CT — iV -)- 1 > fi a + N holds, the difference 
A L (N) = K% ax (\, fi) - n™ ax (\', fjf) is 1 for L = + N, y, a + N + 1, . . . , X a - N+ 1, and zero 
otherwise (with n™ ax as defined in Eq. (0)). Therefore, when A^(A^)=1, a linear combination 
\<p L (N)) = J2 K c n\NHuj(X a - N, ii a + N)kLM) exists such that B \<p L (N)) = 0, and thus our 
claim for family (a) holds. 

The proof for family (b) can be carried out analogously. Here the special irrep (A,/i) = 
(A CT + N, Ha) is obtained by adding N quanta to the 1-direction of the starting configuration. 
In this case there is no restriction on N, hence family (b) is infinite. Note that adding 
quanta to the 3-direction does not yield another family of pure states, since the multiplicity 
for a given L-value in the associated 'special' irreps, (A, /i) = (\ a ,u a — N), decreases as iV 
increases, i.e. Ai(N) < for all L and N. 
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IV. SOLVABLE STATES AND THEIR PROPERTIES 

All OHu states are eigenstates of H(f3 , f3 2 )- They are unmixed and span the entire (X a , fi a ) 
irrep. In contrast, for the excited levels (N > 0), the pure states span only part of the 
corresponding SU(3) irrep. There are other states at each excited level which do not preserve 
the SU(3) symmetry and therefore contain a mixture of SU(3) irreps. 

To construct the pure states for iV > 0, we proceed as follows: Let (A, fi) at Nftw be the 
irrep which contains a pure state with angular momentum L and projection M, \4>lm(N)). 
This state can be written as: 

\4>lm{N)) = c k (L)\NHlu(\, ii)kLM) , (26) 

K=l 

where K™ ax (A,yu) denotes the maximum multiplicity of L in (A,/i), Eq. (0). Obviously, 
\4>lm(N)) is an unmixed eigenstate of H(P ,j3 2 ) if {ip{N — 2)\B \cj) LM (N)) = holds for 
all states \ip(N — 2)) at the (N — 2)ftw level. From the proof it follows that B Q acting 
on states in the 'special' irrep (A, fi) at Nftw can only produce states belonging to the 
'special' irrep (A',//) at (N - 2)Huu, hence ((N - 2)hiu(X', h')k' LM\B \(j) L (N)) = for k' = 
1, 2, . . . , n™ ax (\' i //) ensures that \4>lm(N)) is pure. The ^^(A, /i) coefficients c K (L), which 
characterize the pure state, are thus uniquely determined by the n™ ax (\' , fi') = K™ ax ( A, fi) — 1 
equations 

J2 c n{L) (( - 2) fvuj ( A', fi')K,'LM \B \ Nfvjj{ X, h)kLM) =0 (27) 

K 

and the normalization requirement Y1 K l c «(-^)| 2 = The proof given in the previous section 
guarantees the existence of a solution. 

Making use of the Wigner-Eckart theorem for SU(3) (see Appendix A), the relations in 
Eq. (|27|) can be rewritten as ((A', /x') 1 1 1^ (02) 1 1 1 (A, //)) ^ k c k (L)((A,/x)kL; (02)0|| (A', //)«'£) = 
0, where (_;-||-) denotes a reduced Wigner coupling coefficient for SU(3). Since the triple- 
reduced matrix element of B^ is generally nonzero, we obtain the following conditions: 

c K (L)((\', t i')K'L;(20)0\\(\^)nL) = («' = 1, . . . , ^(A, {,) - 1) . (28) 



K=l 



Note that the matrix elements of the symplectic generators are not relevant for the deter- 
mination of the c K (L), and the SU(3) Wigner coefficients, (_; _||_), can be evaluated numeri- 
cally j|B| or analytically |^ . 



To illustrate the procedure outlined above, we consider the case of 12 C. The leading 
irrep for the nucleus is (A CT , /v) = (0,4), thus the pure states belong to (A, /i) = (0,4) at Ofrw, 
(A,/i) = (2,4) at 2hu, (A, /x) = (4,4) at 4faj, etc. At OJtw, all states (L=0,2,4) are unmixed. At 
2Hu, the possible L-values are 0, 2 2 , 3, 4 2 , 5, 6, and we have Al=o(2^o;)=0 and Al(2Jtw)=1 
for L=2,3,4,5,6. Since the values L=3,5,6 occur only once {KJ£ ax [{2, 4)]=1), the associated 
states are pure (ci(L)=1.0). For L=2,4, for which K™ ax [(2, 4)]=2, the appropriate coefficients 
c K (L) may be determined from the requirements: 

Cl (L)((0,4)L;(2,0)0||(2,4)lL) + c 2 (L)((0, 4)L; (2, 0)0| |(2, 4)2L) = , 

| Cl (L)| 2 + |c 2 (L)| 2 = l . (29) 
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For L=2, one finds ((0,4)2; (2,0)0||(2,4)«2) = -0.85280 (0.05372) for « = 1 (k = 2) g|, and 
thus |0 2M (2^)) = 0.063 |2^(2,4)12M)+ 0.998 \2Hw{2, 4)22M). Similarly, for L=4, one 
obtains ((0,4)4; (2,0)0||(2,4)«4) = -0.75107 (0.23440) for « = 1 (k = 2) [||, and therefore 
Ci(4) = 0.298 and c 2 (4) = 0.955. Analogously, one can proceed for the Atko level. There 
are, for instance, three L=4 states, one of which is pure. One finds: |04m(47jo;)) = -0.637 
|4fiw(4,4)14M)+ 0.761 |4&j(4, 4)24M)- 0.124 |2^j(4, 4)34M), and similarly for the other 
states. 

For a nucleus with (A CT , Ha) = (A, 0), A > 2, pure states with (A',//) = (A — 2,2), 
L = 2, 3, . . . , A — 1, exist at 2huj according to the proof given in Section [ITJ. The odd- 
angular momentum values, L — 3, 5, . . . , A — 1, occur only once (k = 1) and the associated 
states are pure. The even-L values, on the other hand, occur twice, with k = 1 or 2, 
corresponding to Vergados labels and 2, respectively. Since ((A, 0)L; (2, 0)0||(A — 2, 2)kL) 
= [2(A + l) 2 - L(L + l)] 1/2 /[3A(A + 1)] 1/2 for k = 1 and for k = 2 ||, it follows that 
c K (L) = (1.0) for k — 1 (/c = 2). Consequently, the pure K=2 band at 2/jo; consists of 
states with (A',//) = (A - 2, 2), « = 1 (2) for odd (even) L values, i.e. « = n™ ax (X - 2,2). 
An example for such a nucleus is given in Section |VB| , where the 20 Ne system is discussed. 

Having constructed the solvable eigenstates of the PDS Hamiltonian H((3q, /3 2 ), Eq. (p4|), 
we can now give analytic expressions for their energies. We have E(N = 0) = for the Ohw 
level, and 

N 3 
E(N) = f3 2 - (N a - A ff + fi a - 6 + -N) (X a > fi a ) 

N 3 
E(N) = fa-{N a + 2X a + fi a - 3 + -AO (A(j < ^) (30) 

for N > 0. For instance, for N a (X a , fi a ) = 24.5 (0,4), which corresponds to 12 C, this yields: 
E(N = 0) = 0, E{2hu) = 19/3 2 , £(4faj) = 42/3 2 , etc. 

The partial SU(3) symmetry of H(fi ,{3 2 ) is converted into partial dynamical SU(3) sym- 
metry by adding to the Hamiltonian SO(3) rotation terms which lead to L(L + l)-type 
splitting but do not affect the wave functions. The solvable states then form rotational 
bands and since their wave functions are known, one can evaluate the quadrupole transition 
rates between them: 

Z\ 2 5 \{L f \\Q 2 \\U 



B(E2,Lr+L f ) = e 2 b A ( ^ ] - '^"^ Z|ri/I . (31) 

Here b = \Jh/muj is the harmonic oscillator length parameter, Z and A are the nu- 
clear charge and mass, respectively, and the convention for the reduced matrix ele- 
ments is summarized in Appendix A. For unmixed initial and final states, \</>Li(Ni)) = 
J2 Ki c Ki (Li)\Niho(Xi, n^KiLi) and \(j) Lf (N f )) = J2 Kf c Kf {Lf)\Nftko(Xf,fJ,f)K f L f }, the matrix 
element of Q 2 is given by: 

8N t ,Nf s (\i,m)(*f ,n f ) (~ !) ^ V / 6(Csr/3)[(Ai,^)] 

CKi{Li)c Kf (Lf){(\j, HijKiLj] (ll)2||(Ai,/i i )/t / L / ) p=1 

+5 Ni , (Nf+2) V3((\f, H ) 1 1 \A m 1 1 1 (A,, a**)> (32) 
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^c Ki (L,)c K/ (L / )((A i ,/i i )^L i ; (20)2||(A / ,/i / )/« / L / ) 

+^ ll( iv / - 2 )V / 3((A /)/ i / )||| J B( 02 )|||(A,„ / i J )) 
^2c Ki (Li)c Kf (L f )((Xi, ^KiLf, (02)2||(A / ,/i / )K / L / ) , 

where M = for fi = and 1 otherwise. 

For intraband transitions, the above expression reduces to the first term on the right-hand 
side. For interband transitions there are three possibilities: For transitions from Nfiu to 
(N + 2)fkv, the second term has to be evaluated; for Nhu — ► (N — 2)huj transitions, the third 
term is required. For X a ^ 0, /v ^ 0, i.e. for triaxially deformed nuclei, aA r = 0->A r = 
transition is possible as well; in that case the relevant contribution originates from the first 
term. For example, for a transition from Lj = 2 to Lf = in the ground band of 12 C, 
b = 1.668/m, 6<CW 3 )[(0,4)] = 112, and thus B(E2,0frw Li=2^0hiu L f =0) = 0.1925 e 2 /m 4 
x 112/5 x | ((0, 4)2; (1,1)2|| (0,4)0) | 2 = 4.31 e 2 fm 4 = 2.64 W.u. (which corresponds to 4.65 
W.u., when an effective charge e* = 1.327 is used). 

V. APPLICATIONS TO LIGHT NUCLEI 

To illustrate that the PDS Hamiltonians of Eq. (^) are physically relevant, we present 
applications to prolate, oblate, and triaxially deformed nuclei. We compare energy spectra, 
reduced quadrupole transition rates, and eigenstates of 

Hp DS = h(N) + ^H((3 = 12, (3 2 = 18) + l2 L 2 + 74 L 4 (33) 

to those of the symplectic Hamiltonian 

H SpG = Ho- XQ2 ■ Q 2 + d 2 L 2 + d 4 L 4 . (34) 

Here the function h(N), which contains the harmonic oscillator term Hq, is simply a constant 
for a given Nhw excitation. We select light, p-shell and ds-shell, nuclei for which a full, three- 
dimensional symplectic calculation can be carried out, that is, a limitation to a submodel 
of the Sp(6,R) model is not required. Since we employ Hamiltonians composed solely of 
Sp(6,R) generators, we restrict the model space to one Sp(6,R) irrep (represented by one 
'cone' in Fig. 0). We include excitations up to 8fkJ. 

A. The 12 C case 

The first nucleus to be considered is 12 C, with four protons and four neutrons in the 
valence p-shell. This nucleus has previously been studied in the Sp(2,R) submodel of the 
SSM f3l| , [48~1 . Here we employ the full, three-dimensional, symplectic shell model. The 
leading Sp(6,R) irrep for this case is N a (\ ay l^a) = 24.5(0,4). At the 2t\w level SU(3) irreps 
(A,/i) = (2,4), (1,3), (0,2) occur, at the 4hw level we have (A,/i) = (0,6), (1,4), (2,2) 2 , 
(4,4), (3,3), (1,1), (0,0), and so on for higher excitations. The parameters of Hs p g were 
fitted (simultaneously) to the ground band energies and the 2\ — > 0^ reduced quadrupole 
transition strength, for symplectic model spaces including excitations up to 2Jtw, 4fkj, 6hu, 
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and 8ftw. The resulting B(E2) strengths are listed in Table |and several low-lying rotational 
bands are shown in Fig. [3]. The left part of the figure shows the experimental energies of the 
ground band [fj9| , while the center portion (labeled Q 2 ■ Q2) shows the calculated ground 
band (K=0i), as well as several resonance bands which are dominated by 2fko excitations 
(K=2i, O2, li,0s), Ahuj excitations (K=4i), and 6hu> excitations (K=6i). The parameters of 
Hpos were determined as follows: 72 and 74 were fixed by the level splittings of the ground 
band, £ was chosen to fit the energy difference between the K=2i and K=02 bandheads of 
the symplectic calculation, and h(N) was adjusted to reproduce approximately the relative 
positions of the K=0i,2 1 ,4 1 , and 61 bandheads. The resulting spectrum is that shown on 
the right side of Fig. ||, labeled PDS. 

Since Hpp>s does not mix states with different Nftw excitations, the B(E2) values ob- 
tained in the PDS calculations require an effective charge e* = 1.33 to match the experimen- 
tal values (compare Table |). Overall, we find little deviation between the energies 
and electromagnetic transition strengths of the two approaches. A better measure for the 
level of agreement between the PDS and symplectic results is given by a comparison of the 
eigenstates. According to the proof given in Section |TJ, the Hamiltonian Hpp>s should have 
sets of solvable, pure-SU(3) eigenstates, which can be organized into rotational bands: All 
Ohw states should be pure (X a , ji a ) = (0, 4) states, and at 2twj a rotational band with good 
SU(3) symmetry (A, /i) = (2,4) and L = 2,3,4,5,6 should exist. Similarly, we expect pure- 
SU(3) bands at ifau with (A, /i) = (4, 4) and L = 4,5, 6, 7, 8, at 6hw with (A, /i) = (6, 4) and 
L = 6,7, 8, 9, 10, etc. An analysis of the PDS eigenstates shows that this is indeed the case. 
The associated rotational bands are indicated in Fig. |3|. 

Figure f| shows the decomposition of representative (L 71 = 2 + ) states of the five lowest 
rotational bands for the Hs p q and Hpp>s Hamiltonians. The left side of the figure illustrates 
the amount of mixing in the wave functions of the 8fru) (Q2 ■ Q2) calculation: Members of 
the ground band (K=0i) are nearly pure (~ 90%) Ohio states and the resonance bands have 
strong 2/ju; contributions (> 60%). The K=2i, li, and O3 bands contain admixtures from 
Nhu excited states with N > 2, while the K=02 contains admixtures from both the Ohuj 
space and from higher oscillator shells. The relative strengths of the SU(3) irreps within the 
2hu) space are given as well. We find that each rotational band tends to be dominated by 
one representation, namely (2,4) for the K=2i and K=02 bands, (1,3) for K=l 1; and (0,2) 
for K=03, with the other irreps contributing less than 3%. The right side of Fig. |] shows 
the structure of the PDS eigenstates. Since the Hamiltonian Hpp>s does not mix different 
major oscillator shells, each eigenstate belongs entirely to one Nhuj level of excitation. Here 
the ground band belongs to the Ofiw space, while the four resonance bands are pure 2hw 
configurations. Comparing this with the symplectic results, we observe that the Nfiui level 
to which a particular PDS band belongs also dominates the corresponding symplectic band. 
Furthermore, within this dominant excitation, eigenstates of Hs p § and Hpp,s have very 
similar SU(3) structure, that is, the relative strengths of the various SU(3) irreps in the 
symplectic states are approximately reproduced in the PDS case. This holds for the K=0i 
and K=2x bands, which are pure in the PDS scheme, as well as for the mixed K=02, li, and 
O3 bands. The above statements are also true for higher Nfrw excitations, as is illustrated 
in Fig. I for the L = 6 states of the N=2 K=2 1; N=A K=4i, and N=6 K=6i bands. Note 
also that, in the symplectic case, admixtures from higher shells in the L=Q wave functions 
originate predominantly from the 'special' irreps (A,/i) = (N,4). 

The 12 C example given above nicely illustrates the concept of a partial dynamical sym- 
metry for a fermionic many-body system. The pure PDS eigenstates form rotational bands 
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which follow the pattern for solvable states of family (b), their energies and E2 transition 
strengths between them can be evaluated analytically according to Eqs. (|30|)-(|32|). 

B. The 20 Ne case 

We now turn to a system with pure PDS eigenstates that follow pattern (a): The 20 Ne 
nucleus, with two valence protons and neutrons each, has previously been described within 



the symplectic model framework [28|, |43|, 51], |52 |. The leading Sp(6,R) irrep for this prolate 



nucleus is N f7 (X r7} = 48.5(8,0). We expect to find solvable, pure-SU(3) eigenstates of 
Hpos at Ohu), 2ftw, and Ahu. More specifically, there should be a K=0i L = 0,2,4,6,8 
rotational band with (A,/i) = (8,0) at 0ftw, a K=2i L = 2,3,4,5,6,7 band with (A, //) = 
(6,2) at 2hw, and a K=4 1 L = 4, 5 'band' with (A,//) = (4,4) at 4Hu. Pure PDS states at 
higher levels of excitation do not exist. 

As in the 12 C case, we compare the eigenstates of H PDS to those of the symplectic 
Hamiltonian Hs p e- Least squares fits to measured energies and B(E2) values of the ground 
band of 20 Ne |53[ were carried out for 2ftiu, 4hu, 6hu, and 8fiw symplectic model spaces. 
The resulting energies and transition rates converge to values which agree with the data, as 
is illustrated in Fig. ^| and Table 0]. The parameters 72 and 74 in Hpds were determined 
by the energy splitting between states of the ground band, £ was adjusted to reproduce 
the relative positions of the 2hw resonance bandheads and h(N) was fixed by the energy 
difference [.E(O^)— £"(0^ )]. Figure || and Table |J demonstrate the level of agreement between 
the PDS and symplectic results. 

An analysis of the structure of the ground and resonance bands reveals the amount of 
mixing in the 8hu symplectic (Q2 • Q2) wave functions. Figure [7| shows the decomposition 
for representative (L n = 2 + ) states of the five lowest rotational bands. Ground band (K=0i) 
states are found to have a strong Ofru component (> 64%), and three of the four resonance 
bands are clearly dominated (> 60%) by 2hu configurations. States of the first resonance 
band (K=02), however, contain significant contributions from all but the highest Nhw exci- 
tations. The relative strengths of the SU(3) irreps within the 2hu space are shown as well: 
as in the 12 C case, states are found to be dominated by one representation [(10,0) for the 
K=0 2 band, (8,1) for K=li, (6,2)k = 2 for K=2 U and (6,2)k = 1 for K=0 3 here], while the 
other irreps contribute only a few percent. Such trends are present also in the more realistic 



symplectic calculations of Suzuki [52 



As expected, Hpds has families of pure-SU(3) eigenstates which can be organized into 
rotational bands, Fig. []. The ground band belongs entirely to N = 0, (A,/i) = (8,0), and 
all states of the K=2 X band have quantum labels N = 2, (A,/i) = (6,2), k = 2. The K=4x 
band at Ahw is not shown. A comparison with the symplectic case shows that the Nfkv level 
to which a particular PDS band belongs is also dominant in the corresponding symplectic 
band, Fig. [7|. As before, within this dominant excitation, eigenstates of Hpds and Hs p e 
have similar SU(3) distributions; in particular, both Hamiltonians favor the same (A,/i)/t 
values. Significant differences in the structure of the wave functions appear, however, for the 
K=02 resonance band. In the 8hu symplectic calculation, this band contains almost equal 
contributions from the Ohu, 2hw, and 4hw levels, with additional admixtures of 6Hlj and 
8^0; configurations, while in the PDS calculation, it belongs entirely to the 2hu level. These 
structural differences are also evident in the interband transition rates, as is illustrated 
in Table |TJ Whereas the intraband B(E2) strengths are approximately equal in both 
approaches, we observe that the interband rates differ by a factor of 2-3 in most cases. 
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These differences reflect the action of the inter-shell coupling terms that are present in the 
quadrupole-quadrupole interaction of Eq. (p3|), but do not occur in the PDS Hamiltonian. 
Increasing the strength \ of Q2 ■ Qi in H$ p % will also spread the other resonance bands over 
many Nftw excitations. The K=2i band (which is pure in the PDS scheme) is found to resist 
this spreading more strongly than the other resonances. For physically relevant values of x, 
the low-lying bands have the structure shown in Fig. [7]. 



C. The 24 Mg case 

The final example to be considered here involves the triaxially deformed nucleus 24 Mg, 



which has been the subject of several symplectic model studies [|42], [43], |56|, p7[| . With four 
valence protons and neutrons in the ds-shell each, and a leading Sp(6,R) irrep N a (\ a ,fi a ) 
= 62.5(8,4), this system is the most complicated one to be investigated here. Since both 
\ a 7^ and \i a 7^ 0, the symplectic Hilbert space has a very rich structure. The (8,4) 
representation at Ohw contains three rotational bands: a K=0 band with L — 0,2, 4, 6, 8, a 
K=2 band with L = 2, 3, . . . , 10, and a K=4 band with L = 4, 5, . . . , 12. At the 2hu level, 
there are six possible SU(3) irreps, (10,4), (8,5), (6,6), (9,3), (7,4), and (8,2), which contain a 
total of four K=0, two K=l, four K=2, two K=3, three K=4, one K=5, one K=6, . . . bands. 
At the 4hu level, there are 15 different SU(3) irreps, at 6hu>, there are 25, etc. Accordingly, 
the number of states for a given angular momentum value L increases dramatically with the 



inclusion of higher excitations. This is illustrated in Table [IV| for L — 0,1,2, ... ,8. 

Since the interactions in Hs p e do not distinguish different ^-multiplicities, it becomes 
necessary to make use of the integrity basis operators A3 and A 4 discussed in Section II, 
which allow us to reproduce the experimentally observed K-band splitting in the spectrum 
of 24 Mg. Using the Hamiltonian 

H' Sp6 = H Sp6 + C3A3 + C4A4 , (35) 

we carried out least squares fits to measured energies and B(E2) values for 2hw, 4ftw, and 
6huu symplectic spaces. 

Figure |8| (top) displays the energies obtained with the 6hj calculation (right part of 
the figure) in comparison with the experimental values J5B| , |5T| (left side). In addition to 
the ground (K=0i) and 7 (K=2i) bands, the calculated K=4i band, which is dominated 
by Oftw configurations, and several low-lying symplectic K=0, 2, 4, and 6 bands, which are 
dominantly 2frw resonances, are shown. Table V lists various B(E2) transition rates between 
the low-lying states of 24 Mg. We find that the results of the symplectic calculations are in 
good agreement with the data 4 . Specifically, the 7-band is correctly located and nearly all 
the calculated intraband and interband transition rates fall, without the use of an effective 
charge, within experimental uncertainties. The Ahu results are better than the 2hu results, 
with the 6hu calculation yielding only moderate improvements. 



4 Note that we have used the experimental B(E2) values from Ref. |5£|, which contains a more complete 
list of B(E2) data than the compilation by Endt |39|. The latter gives values of 20.6±0.4 W.u., 35±5 
W.u., and 37±12 W.u. for the first three transitions listed in Table V. Fitting the symplectic Hamiltonian 
parameters to reproduce the values of Ref. |5{| gives results very similar to the ones presented here and 
does not alter our conclusions. 
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In analogy with the symplectic case, we include terms A3 and X4 in the PDS Hamiltonian: 

H'pds = H PDS + C3X3 + c 4 X 4 . (36) 

As we will see below, the introduction of these extra terms breaks the partial symmetry. We 
fixed C3 and C4 at the values that were used in the 6hw symplectic calculation, determined 72 
and 74 from the level splittings in the K=0i ground band, and adjusted £ so as to reproduce 
the relative positions of selected 2hu bandhead states (we focused on the lowest three K=0 
bands and the first K=6 band). Then h(N) was chosen to reproduce approximately the 
positions of the 2hu resonances relative to the ground and 7 bands. 

We obtain an energy spectrum which agrees well with the results of the symplectic cal- 
culation, as is shown in Fig. ||. The B(E2) strengths for the ground and 7-bands, rescaled 
by an effective charge e*=1.75, are given in Table V. We find good agreement between 
the PDS and symplectic calculations for the intraband transitions, whereas there are larger 
deviations in the interband rates. 



According to the proof given in Section [II, the three rotational bands at Ohio should be 



pure in the PDS scheme, and at 2hw there should be a (short) rotational K=6 band with 
L—6,7, which belongs entirely to the (A,/i) = (6,6) representation. We find that the Ohu 
states are indeed pure, but the K=6 L = 6, 7 band has small admixtures (< 1%) from 2huj 
irreps other than (A,/i) = (6,6), thus indicating that H' PDS , unlike Hpas, is not an exact 
partial dynamical symmetry Hamiltonian, due to the presence of the K-band splitting terms 
X3 and X4. This can be understood as follows: While X3 and X4 cannot mix different 
SU(3) irreps, their eigenstates involve particular linear combinations of different k values. 
Since the operators X3 and X4 do not commute with Bq, these linear combinations differ 
from configurations resulting from the PDS requirement Bq\4>) = 0. Fortunately, a very 
small amount of symmetry-breaking suffices to fit the relative positions of the ground and 
7-bands, as can be inferred from the eigenstate decompositions plotted in Fig. [| Shown 
are the decompostions of the L = 6 states associated with the calculated H' Sp6 and H' PDS 
spectra. More specifically, we have plotted the contributions from the SU(3) irreps at Ohu 
and 2hw, as well as the (summed) contributions from 4hu and Qhw excitations. 

As in the previous examples, we observe that the eigenstates of both Hamiltonians have 
very similar structures: For a given state, the same Nhuj level of excitation is dominant in 
both calculations and, moreover, within this dominant excitation, we find similar SU(3) dis- 
tributions. The structural differences that do exist are, again, reflected in the very sensitive 



interband transition rates, as can be seen in Table VI 



VI. COMPARISON OF PARTIAL SYMMETRIES IN BOSONIC AND 
FERMIONIC MANY-BODY SYSTEMS 

Partial dynamical symmetries were first studied in the Interacting Boson Model (IBM) of 
nuclei . In [[16] , the following IBM Hamiltonian was used to reproduce measured energies 
and E2 rates of 168 Er: 

H IBM (h , h 2 ) = h P£Po + h 2 P\ ■ P 2 , (37) 
where ho, h 2 are arbitrary parameters and P L , L = and 2, are boson pair operators: 

pt = S ■ S - 2(s f ) 2 , 

Pi = 2sUl + >/7(dW)« . (38) 
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The creation operators s' and (// = 0, ±1,±2) denote a monopole boson with angular 
momentum and parity J n = + , and a quadrupole boson with J 71 " = 2 + , respectively. They 
represent correlated valence nucleon pairs and are the basic building blocks of the IBM. 
The pair operators Pq and PL are components of a (A,/i) = (0,2) SU(3) tensor, and their 

Hermitean adjoints, P and P 2ji = (— l) M p2 ) _ At , are characterized by (A,/i) = (2,0). 

It can be shown that for h 2 = ho, the Hamiltonian of Eq. fl37|) becomes a SU(3) scalar 
(related to the Casimir operator of SU(3)) and for h 2 = —ho/5, it transforms as a (A,/i) = 
(2,2) SU(3) tensor component. In general, HiBM^ho, h 2 ) is therefore not a SU(3) scalar, 
nevertheless it turns out that it always has an exact zero-energy eigenstate, denoted in what 
follows by |c; N), where the integer N gives the total number of bosons in the system. The 
state \c; N) describes a condensate of bosons and can be written as 

1 r i N 
|c; N) = -j= + V24)/V3\ |0) . (39) 

It is the lowest weight state in the SU(3) irrep (A,/i) = (2iV, 0) and serves as an intrinsic 
state for the SU(3) ground band. The rotational members of the ground band with good 
angular momentum L are obtained by projection from \c;N). Moreover, one finds that 
states of the form 

\k)^{Pl 2 ) k \c-N) (40) 

are eigenstates of H IBM (h , h 2 ) with eigenvalues E k = Qh 2 {2N + 1 — 2k)k and good SU(3) 
symmetry (2N — 4k, 2k), where 2k < N. They are lowest weight states in these representa- 
tions and serve as intrinsic states representing 7 fc bands with angular momentum projection 
K = 2k along the symmetry axis. 

Since HjBM{ho, h 2 ) is rotationally invariant, it follows that states of good L projected from 
\k = 0) = \c;N) and \k), k > 0, are also eigenstates with energy Ek and SU(3) symmetry 
(2N — 4k, 2k). The projected states span the entire (2N,0) representation for k = 0, but 
only part of the corresponding irrep for k > 0. There are other excited states which do not 
preserve the SU(3) symmetry and therefore contain a mixture of SU(3) irreps, including the 
'special' irreps {2N—4k, 2k). Since H IBM (h , h 2 ) is not a SU(3) scalar, but possesses a subset 
of solvable eigenstates with good SU(3) symmetry, it is a partial symmetry Hamiltonian. 
Adding L 2 , the Casimir operator of SO(3), to H IBM (h , h 2 ) converts the partial symmetry 
to a partial dynamical symmetry and contributes a L{L + 1) splitting, but does not affect 
the wave functions. 



The boson and fermion Hamiltonians, HjBM(ho, h 2 ) of Eq. (|37|) and H((3o, j3 2 ) of Eq. (p~L|) . 
have several features in common: Both display partial SU(3) symmetry, they are constructed 
to be rotationally invariant functions of (A,/i) = (2,0) and (A,/i) = (0,2) SU(3) tensor 
operators, and SU(3) tensor decompositions show that both contain (A, /i) = (0, 0) and (2,2) 
terms only. H IBM (h ,h 2 ), as well as H(j3 ,j3 2 ), has solvable pure-SU(3) eigenstates, which 
can be organized into rotational bands; the degeneracies within these bands are lifted by 
adding the SO (3) term L 2 to the Hamiltonian. The ground bands are pure in both cases, 
and higher-energy pure bands coexist with mixed-symmetry states. 

There are several significant differences between the bosonic and fermionic PDS Hamilto- 
nians, however. For example, the ground band of the Hamiltonian Hi B M(ho, h 2 ), Eq. (0), is 
characterized by (A, /i) = (2N, 0), i.e., it describes an axially-symmetric prolate nucleus. Is is 
also possible to find an IBM Hamiltonian with partial SU(3) symmetry for an oblate nucleus. 
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It can be shown that these two cases exhaust all possibilities for partial SU(3) symmetry 
with a two-body Hamiltonian in the IBM-1 with one type of monopole and quadrupole 
bosons. In contrast, the fermionic Hamiltonians considered here can accomodate ground 
bands of prolate [(A CT ,0)], oblate [(0, //„■)], and triaxial [(X a ,/i a ) with A CT 7^ 0, /i a 7^ 0] shapes. 

Another difference between the fermionic and the bosonic PDS cases discussed here lies 
in the physical interpretation of the excited solvable bands. While these bands represent 7, 
double-7, etc. excitations in the IBM, they correspond to giant monopole and quadrupole 
resonances in the fermion case. 

Furthermore, whereas the pure eigenstates of H IBM (h , h 2 ) can be generated by repeated 
action of the boson pair operator P\ 2 on the boson condensate and subsequent angular mo- 
mentum projection, a similar straightforward construction process for the special eigenstates 
of H(f3 ,j3 2 ) has not been identified yet. The situation seems to be more complicated in the 
fermion case, which is also reflected in the fact that H(f3 ,f3 2 ) has two possible families of 
pure eigenstates, one finite, the other infinite. The association of the special states to one 
or the other family depends on the OHuj symplectic starting configuration. 

The comparison of partial dynamical symmetries in bosonic and fermionic systems above 
illustrates that, in spite of similar algebraic structures of the associated Hamiltonians, two 
given systems with partial symmetries may exhibit not only different physical interpreta- 
tions, but also different systematic features and different mechanisms for generating the 
partial symmetries in question. 

VII. SUMMARY AND CONCLUSIONS 

The fundamental concept underlying algebraic theories in quantum physics is that of 
an exact or dynamical symmetry. Realistic quantum systems, however, often require the 
associated symmetry to be broken in order to allow for a proper description of some observed 
basic features. Partial dynamical symmetry describes an intermediate situation in which 
some eigenstates exhibit a symmetry which the associated Hamiltonian does not share. The 
objective of this approach is to remove undesired constraints from the theory while preserving 
the useful aspects of a dynamical symmetry, such as solvability, for a subset of eigenstates. 

We have presented an example of a partial dynamical symmetry in an interacting many- 
fermion system. In the framework of the symplectic shell model, we have constructed a family 
of rotationally invariant Hamiltonians with partial SU(3) symmetry. We have demonstrated 
that the PDS Hamiltonians are closely related to the deformation-inducing quadrupole- 
quadrupole interaction and break SU(3) symmetry, but still possess a subset of 'special' 
solvable eigenstates which respect the symmetry. The construction process for these special 
states was outlined and analytic expressions for their energies and for E2 transition rates 
between them were given. 

To illustrate that the PDS Hamiltonians introduced here are physically relevant, we have 
presented applications to oblate, prolate, and triaxially deformed nuclei. Specifically, we 
have compared the energy spectra, reduced quadrupole transition strengths, and eigenstate 
structures of the partial symmetry Hamiltonians to those of a symplectic shell model Hamil- 
tonian containing a realistic quadrupole-quadrupole interaction. Although the PDS Hamil- 
tonians cannot account for intershell correlations, we have observed that various features of 
the quadrupole-quadrupole interaction are reproduced with a particular parameterization 
of the partial symmetry Hamiltonians. PDS eigenfunctions do not contain admixtures from 
different Nfiw configurations, but belong entirely to one level of excitation. We have found 
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that, for reasonable interaction parameters, the NJtuj level to which a particular PDS band 
belongs is also dominant in the corresponding band of exact Q2 ■ Q2 eigenstates. Moreover, 
within this dominant excitation, eigenstates of both Hamiltonians have similar SU(3) distri- 
butions. Structural differences, nevertheless, do arise and are reflected in the very sensitive 
interband transition rates. Overall, however, we may conclude that PDS eigenstates approx- 
imately reproduce the structure of the exact Q 2 ■ Q2 eigenstates, for both ground and most 
resonance bands. 

The notion of partial dynamical symmetries extends and complements the familiar con- 
cepts of exact and dynamical symmetries. It is applicable when a subset of physical states 
exhibit a symmetry which does not arise from the invariance properties of the relevant 
Hamiltonian. Recent studies, including the one presented here, show that partial symme- 
tries may indeed be realized in various quantum systems. This indicates that PDS is not 
a mere mathematical concept, but may serve as a practical tool in realistic applications of 
algebraic methods to physical systems. 
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Appendix A: SU(3) Wigner coefficients and Wigner-Eckart theorem 



If a represents a set of labels used to distinguish orthonormal basis states 
within a given irreducible SU(3) representation (A,/i), the Wigner coefficients 
( (Aj, Hi)oli, (A 2 , ^2)^2 I (A, n)(x)p are defined as the elements of a unitary transformation 
between coupled and uncoupled orthonormal irreps of SU(3) in the a-scheme |37fl: 



\(X,ii)a) p = (Ai,/xi)ai; (A 2 ,/i 2 )a 2 | (\,p)a) p | (Ai,/ii)ai )| (A 2 ,/i 2 )a 2 ), (41) 

OL10L2 

and the inverse transformation is given by: 

I (Ai,/ii)ai)| (A 2 ,/i 2 )a 2 ) = ( (Ai, pxjar, (A 2 , p 2 )a 2 1 (A, p)a ) p |(A, p)a ) p . (42) 

p(\,fj,)a 

Here a = eAM\ for the SU(3) D SU(2) <g) U(l) (canonical) group chain and a = nlm for 
the SU(3) D SO(3) reduction employed in this work. The subgroup chains impose certain 
restrictions on the above couplings, for example the usual angular momentum coupling rules, 
I = h + h, ■ ■ ■, \h — h\, and m = mi + m 2 apply for the chain containing SO (3). 

The outer multiplicity label p = 1,2,..., p max is used to distinguish multiple occurrences 
of a given (A, p) in the direct product (Ai,/ii) x (A 2)/ u 2 ): p = 1, 2, . . . , p max , where p max 
denotes the number of possible couplings (A 1; //i) x (A 2 ,/x 2 ), and the possible (A, p) irreps in 
the product can be obtained by coupling the appropriate Young diagrams [||]. O'Reilly |60|| 
determines a closed formula for the decomposition of the outer product (Xi,p\) x (A 2 ,/x 2 ) 
of SU(3) irreps for arbitrary positive integers Aj, pi and derives necessary and sufficient 
conditions for a SU(3) irrep (A,/i) to appear as summand in (Ai,/xi) x (A 2 ,/x 2 ). 



23 



It is possible to factor out the dependence of the above SU(3) D SO(3) Wigner coupling 
coefficient on the m subgroup label by defining so-called double-barred or "reduced" SU(3) 
coupling coefficients: 

( {Xi^i) Kxhmi, (A 2 ,// 2 )/« 2 Z 2 m 2 1 (A,//)«Zm) p 

= ( (Ai,/xi)ki/i; (A 2 ,/i 2 )/t 2 / 2 || (\,h)kI) p (km 1 , l 2 m 2 \lm) . (43) 

s v ' v * ' 

reduced Wigner coefficient geometric part 

The "geometric" part (__|_) is simply a SU(2) Clebsch-Gordan coefficient. From the unitarity 
of the full SU(3) Wigner and the ordinary SU(2) Clebsch-Gordan coefficients it follows 
that the double-bar coefficients are also unitary. With the phase convention introduced 
in Ref. [57J they become real, and therefore orthogonal. Draayer and Akiyama j?7| give a 



prescription for the unique determination, including the phases, of SU(3) Wigner coefficients 
and derive their relevant conjugation and symmetry properties. They furthermore provide a 
computer code which allows for a numerical determination of the coefficients |^6[|. Analytic 



expressions for Wigner coefficients which are of particular interest in p-shell and ds-shell 
nuclear shell-model calculations are tabulated in Ref. |JT) for the canonical subgroup chain 
and in Ref. |3]| g7|, §2| for the SU(3) D SO (3) chain. 



The Wigner-Eckart theorem for the group SU(2) yields SU(2)-reduced (double-bar) ma- 
trix elements of a SO (3) irreducible tensor operator: 

(l 3 m 3 \T l ^\l imi ) 

= ( km i; l 2 m 2 \ l 3 m 3 ) ^4f=^ ■ ( 44 ) 

Analogously, the generalized Wigner-Eckart theorem allows one to express matrix elements 
of SU(3) irreducible tensor operators as a sum over p of the product of a p-dependent 
generalized reduced matrix element multiplied by the corresponding Wigner coefficient |57|| : 



((X 3 ,p 3 )a 3 \T^^\(X 1 , f i 1 )a 1 ) 

= J2 ( CWi) «i; ( A 2, // 2 ) a 2 \ (A 3 , /i 3 )« 3 )p ( (A 3 , /is) 1 1 1 T^) 1 1 |(Ai, /ii) ) p . (45) 
p 

For more details on SU(3) coupling and recoupling coefficients, see the compilation in Ap- 
pendix C of Ref. |34| and references therein. 



Appendix B: Matrix elements of relevant operators 

The calculations presented here require expressions for matrix elements of the Sp(6,R) 
generators A^ 20 \ B^ 02 \ and C^ ll \ and combinations thereof. None of these operators con- 
nect states belonging to different symplectic representations and, furthermore, the SU(3) 
generators = L q and C 2 ^ = -^Q 2fl act only within one level of excitation, N. Matrix 
elements for C^ 11 -* in the standard SU(3) bases are given by [[H], 61 



((AV)|||£ (11) |||(A,/i)> = (-l)*^2(Csu3)[(\,ri)5 { w M , (46) 
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where Csu3 denotes the second-order Casimir operator of SU(3), given in Eq. (|), and 
M = 1 for fi ^ and ^ = for fi = 0. The reduced matrix element ((X' , /j')\\\C^\\\(X, jj)) 
is related to the full SU(3) matrix element via the Wigner-Eckart theorem for SU(3) and 
the phase is chosen to be consistent with that of Ref. [44|. 



Several strategies for calculating matrix elements of the symplectic generators and 
have been explored. A direct way is to use the Sp(6,R) commutation relations to 
derive recursion formulae, as shown by Rosensteel []33]. Another approach is to start from 
approximate matrix elements and to proceed by successive approximations, adjusting the 
matrix elements until the commutation relations are precisely satisfied PBfl . Deenen and 
Quesne p3 have employed a boson mapping to obtain generator matrix elements, and 
Castanos et al. [|64j have derived simple analytical functions for some special irreps. The 



most elegant method, outlined by Rowe in Ref. p5| , involves vector- valued coherent state 
representation theory and evaluates matrix elements of the symplectic raising and lowering 
operators by relating them to the matrix elements of a much simpler u(3) ® Weyl algebra. 
A listing of the relevant formulae is beyond the scope of this appendix, the reader is thus 



referred to Ref. |65[ for details of the calculation. 

Matrix elements of the SU(3) D SO(3) integrity basis operators X 3 = (L x Q E )(\) ■ L and 
X± = (L X Q E )(i) ■ {L x Q E )(i) can be given in terms of SO (3) Racah recoupling coefficients 
Z 2 , Ua; W23) H and the SU(3) D SO(3) reduced matrix elements of 0: 



((A, ^)Klrn\X 3 \(X' , fi')n'l'm') 

= 0~(X', f i>)(X, f i)0~l'lO~m'm3l(l + l)V / 2l + 1 

xW(l, 1, 1, 1; 1, 2)((A, h)k1\ \C^\ |(A, ^k'I); (47) 
((A, fi)Klm\X4\(X' , n')K'l'm') 

= ^(X'^'XX^l'lSm'm^K 1 + l)v2Z + 1 

x ^(-l) w " +1 v / 27 7 TT [W(l, 1, 2, 1"; I, l)] 2 

k"1" 

x ((A, p)Kt\ \Ci n) \ |(A, ii)k"1")({X, ^)k"1"\ \Ct ] I |(A, ^)k'1). (48) 
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TABLE I: B(E2) values (in Weisskopf units) for ground band transitions in 12 C. Compared are sev- 
eral symplectic calculations, PDS results, and experimental data. Q denotes the static quadrupole 
moment of the L n = 2\ state and is given in units of eb. The experimental values are taken 
from Refs. |5(|. PDS results are rescaled by an effective charge e*=1.33 and the symplectic 
calculations employ bare charges. 

Transition Model B(E2) [W.u.] B(E2) [W.u.] 
Ji Jf 2Hlu AHuj 6Hlu Shw PDS Exp. 
2^ 4.65 4.65 4.65 4.65 4.65 4.65 ± 0.26 
4^ 2 4.35 4.27 4.24 4.23 4.28 n/a 
Q [eb] 0.059 0.060 0.060 0.060 0.058 0.06± 0.03 



TABLE II: B(E2) values (in Weisskopf units) for ground band transitions in 20 Ne. Compared are 
experimental data, predictions from several symplectic calculations, and PDS results. The static 
quadrupole moment of the L n = 2~f state is given in the last row. The experimental values are 



taken from Refs. |53| , |54| , p5f |. PDS transition rates are rescaled by an effective charge e*=1.95, 
while the symplectic calculations use bare charges. 

Transition Model B(E2) [W.u.] B(E2) [W.u.] 



Ji- 




2hw 




6fkv 


8hu) 


PDS 


Exp. 


2 - 


-> 


14.0 


18.7 


19.1 


19.3 


20.3 


20.3 ± 1.0 


4 - 


-> 2 


18.4 


24.5 


24.6 


24.5 


25.7 


22.0 ± 2.0 


6 - 


-> 4 


17.1 


22.3 


21.5 


20.9 


21.8 


20.0 ± 3.0 


8 - 


-> 6 


12.4 


15.2 


13.3 


12.4 


12.9 


9.0 ± 1.3 



Q [eb] -0.14 -0.16 -0.16 -0.16 -0.17 -0.23 ± 0.03 



28 



TABLE III: Comparison of intraband and interband B(E2) rates for 20 Ne. Shown are various 
transitions between states of the lowest rotational bands. K=0i denotes the ground band, which 
is dominated by OHuj configurations; members of the other bands correspond to 2huj resonances. 
Results are from the PDS calculation (rescaled by e*=1.95) and from the 8ftw symplectic approach 
(e* = 1.0). In the last column, ratios of the calculated transition strengths are given. 

Transition Model B(E2) [W.u.] BE2(PDS) 



Ji Ki J f K f 




PDS 




2 Oi 





Oi 


1 Q 3 


9(1 3 


1 


2 Do 
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Di 
u i 






9 1 fi 


2 3 





Oi 


U. 1U 


n ^9 


o. ID 


2 Oi 





o 2 


2.9 


5.7 


1.94 


2 2 





o 2 


20.3 


27.8 


1.37 


2 3 





o 2 


0.15 


0.13 


0.84 


2 Oi 





o 3 


0.17 


0.48 


2.80 


2 2 





o 3 


0.25 


0.26 


1.01 


2 3 





o 3 


12.9 


16.8 


1.30 


4 Oi 


2 


Oi 


24.5 


25.7 


1.05 


4 2 


2 


Oi 


10.9 


22.8 


2.09 


4 li 


2 


Oi 


2.3 


5.8 


2.55 


4 2x 


2 


Oi 


0.63 


2.3 


3.66 


4 3 


2 


Oi 


0.09 


0.30 


3.34 
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TABLE IV: Dimensions of symplectic Hilbert spaces for Mg. Shown are the number of L-states 
(L = 0, 1, . . . , 8) for spaces which include Nftw excitations up to N = 0, 2, 4, and 6, respectively. 

Symplectic Angular momentum L 

space 1 2 3 1 5 6 7 8~ 

mo Tl) 2 1 3 2 3 2 3 

(0+2)^j 4 3 11 10 17 15 19 16 18 
(0+2+4)^; 13 15 40 41 62 59 71 63 67 
(0+2+4+6)/ilj 32 49 110 122 171 171 198 182 187 
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TABLE V: B(E2) strengths of 24 Mg. Compared are results from 2Huj, Ahu), and 6Hu> symplectic 
calculations, a PDS calculation, and experiment [55|, 58]. Both intraband and interband transitions 



between states of the ground (K=0i) and 7 (K=2i) band are given. The static quadrupole moment 
of the 2^~ state is listed in the last line (in units of eb) a . The symplectic model reproduces the 
observed transition rates without employing effective charges, while the PDS approach requires 
e*=1.75. 

Transition Model B(E2) B(E2) 



Ji Ki 






2huj 




6huj 


PDS 


Exp. 


2 


Oi 





Oi 


17.2 


20.2 


20.4 


20.5 


20.5±0.6 


4 


Oi 


2 


Oi 


24.5 


26.9 


26.9 


26.2 


23±4 


6 

8 


Oi 
Oi 


4 

6 


Oi 
Oi 


25.2 
24.4 


25.5 
19.4 


25.2 
19.2 


22.5 
13.6 


34±?g 
16±f 


3 


2i 


2 


2i 


31.6 


35.6 


35.3 


36.6 


34±6 


4 


2i 


2 


2i 


9.7 


11.2 


11.0 


11.6 


16±3 


5 


2i 


3 


2i 


15.3 


17.0 


16.6 


16.8 


28±5 


5 


2i 


4 


2i 


17.3 


18.0 


17.7 


18.0 


14±6 


6 


2i 


4 


2i 


15.3 


19.4 


18.3 


20.1 


23±f 


8 


2i 


6 


2i 


12.4 


18.0 


15.9 


19.6 


>3 


2 


2i 





Oi 


1.1 


1.3 


1.3 


3.1 


1.4±0.3 


2 


2i 


2 


Oi 


2.2 


1.7 


1.9 


3.4 


2.7±0.4 


3 


2i 


2 


Oi 


1.9 


2.4 


2.3 


5.6 


2.1±0.3 


4 


2i 


2 


Oi 


0.2 


1.0 


0.9 


2.7 


1.0±0.2 


4 


2i 


4 


Oi 


2.9 


2.1 


2.3 


4.1 


l.Orbl.O 


5 


2i 


4 


Oi 


1.0 


2.4 


2.0 


6.0 


3.9±0.8 


6 


2i 


4 


Oi 


0.2 


1.2 


1.0 


3.2 


0-8±°;I 




<2 


[eb] 




-0.171 


-0.186 


-0.185 


-0.191 


-0.18±0.02 



a Measurements have given results for |Q| ranging from less than 0.16 eb to nearly double that value. We 
list the value adopted in the review by Spear |55|. 
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TABLE VI: Comparison of intraband and interband B(E2) rates for 24 Mg. Shown are selected 
transitions between states of the K=0i, 02,03, and 22 bands. The PDS values are rescaled by 
e*=1.75. Ratios of the results from the two theoretical approaches are given in the last column. 

Transition Model B(E2) [W.u.] BE2(PDS) 



Ji 




Jj K f Sp(6,R) 
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BE2fSo6) 
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Ol 
u l 
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FIG. 1: Basis construction in the symplectic model. SU(3)-coupled products of the raising 
operator ^4( 20 ) with itself act on an Elliott starting state with (\ a ,^ a ) = (0,[i) ({<7i,<72 = ^1,^3}) 
to generate symplectic 2fou>, Ahw, . . . excitations. Also shown are the SU(3) labels (A, /x) and quanta 
distributions {u\, 0^2,^3} for some excited states. 




(x l hi) a 2 n 2 ) a 3 n 3 > a 4 n 4 ) 



FIG. 2: Symplectic shell model space. The schematic plot illustrates a model space with multiple 
symplectic representations. Each 'cone' corresponds to a Sp(6,R) irrep and is uniquely charac- 
terized by U(3) quantum numbers N a (X a , fj, a ), where (X a ,fi a ) denotes the Elliott SU(3) quantum 
labels for the associated Oftw shell model configuration. For a given starting representation (X a , /j, a ) 
(a = 1,2,3,4 here), one obtains multiple SU(3) configurations, (Xu,,fi w ), at each Nftw level of ex- 
citation (N > 0), indicated here by small filled circles. 
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FIG. 3: Energy spectra for 12 C. Comparison between experimental values (left) [^9|, results from 
a symplectic 8hio calculation (center) and a PDS calculation (right). K=0i indicates the ground 
band in all three parts of the figure. In addition, resonance bands dominated by 2hu> excitations 
(K=2i, O2, li, O3), Afojj excitations (K=4i), and 6Huj excitations (K=6i) are shown for the Sp(6,R) 
and PDS calculations. Additional mixed resonance bands (not shown), dominated by ^frw and 
6^0; excitations, exist for this nucleus. The angular momenta of the positive parity states in the 
rotational bands are L=0,2,4,. . . for K=0 and L=K,K+l,K+2, . . . otherwise. Bands which consist 
of pure-SU(3) eigenstates of the PDS Hamiltonian are indicated. 
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FIG. 4: Decompositions for calculated L 71 " = 2 + states of C. Individual contributions from the 
relevant SU(3) irreps at the Ohuj and 2hw levels are shown for both a symplectic 8fuv calculation 
(denoted Q2 ■ Q2) and a PDS calculation. In addition, the total strengths contributed by the NHuj 
excitations for N > 2 are given for the symplectic case. 
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FIG. 5: Decompositions for calculated L 71 = 6 + states of 12 C. The structures shown are repre- 
sentative for the members of the K=2i,4i, and 61 rotational bands, respectively. States of these 
bands are dominated by NHu excited configurations with (A, fx) = (N, 4), N = 2,4,6,8, in the 
symplectic scheme and are pure in the PDS approach. 
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FIG. 6: Energy spectra for Ne. Experimental ground band (K=0i) energies |53| are shown on 
the left, while theoretical results for both the ground band and 2hu) resonances (K=02, li, 21,03) 
are given in the center and on the right, for a symplectic 8?kv and a PDS calculation, respectively. 
Rotational bands which consist of pure eigenstates of the PDS Hamiltonian are indicated. 
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FIG. 7: Decompositions for calculated L 71 " = 2 + states of 20 Ne. Individual contributions from the 
SU(3) irreps at the Oftw and 2hu> levels are shown for both a symplectic 8hu> calculation (left side) 
and a PDS calculation (right side). For the symplectic approach the summed contributions from 
SU(3) irreps at higher (N > 2) excitations are given as well. 
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FIG. 8: Energy spectra for Mg. Energies from a PDS calculation (bottom) are compared to 
symplectic results (top). Both O/iw-dominated bands (K=0i, 2±, 4i) and some 2ftw resonance bands 
(K=02, 03, 04, 22, 23, 42, 43, 61) are shown. The K=0i,2i,4i (61) states are pure (approximately 
pure) in the PDS scheme. Experimental values for the ground and 7-band energies, taken from 
Refs. |58|, 59 1, are given on the left. 
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FIG. 9: Decompositions for calculated L n = 6 + states of 24 Mg. Eigenstates resulting from the 
symplectic Qhu; calculation are decomposed into their Ofko, 2hw, 4hu>, and 6huj components (denoted 
by Q2 Q2 in the figure). At the Ohio and 2fko levels, contributions from the individual SU(3) irreps 
are shown, for higher excitations (N > 2) only the summed strengths are given. Eigenstates of the 
PDS Hamiltonian belong entirely to one Ntvjj level of excitation, here Ohuj or 2hu>. Contributions 
from the individual SU(3) irreps at these levels are shown. Members of the K=0i, 2i, 4i bands are 
pure in the PDS scheme, and K=6i states are very nearly (> 99%) pure. 



